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1. Introduction 

Fifty years ago Casimir |I| showed that the presence of boundaries changes the ground 
state of the electromagnetic field, leading to nontrivial quantum effects like the Casimir 
force (see also [2-4]). Last years, an attention of many authors [5-40] was attracted to 
the nonstationary modifications of the Casimir effect in the case of moving boundaries 
(a detailed list of publications before 1995 was given in |2l|). The present article is 
devoted to the special case of the nonstationary Casimir effect (NSCE), namely, to the 
effect of photon creation from vacuum in an ideal one-dimensional cavity (a model of 
the Fabry-Perot interferometer) with vibrating boundaries. 

As was understood recently [10,14-23], notwithstanding that the maximal velocity 
of the boundary achievable under the laboratory conditions is very small in comparison 
with the speed of light, a gradual accumulation of the small changes in the quantum 
state of the field could result finally in a significant observable effect, if the boundaries 
of a cavity perform small oscillations at a frequency uo w which is an integer multiple of 
the unperturbed eigenfrequency of the fundamental electromagnetic mode uj\ = irc/L 
(where L is the mean distance between the walls): uj w = pu\, p = 1, 2, . . . (remember 
that the spectrum of the electromagnetic modes is equidistant in the case involved: 
the unperturbed frequency of the pth mode equals u p = pu>i). The time evolution of 
the field in the short time limit suit <C 1 (where e <C 1 is a ratio of the amplitude 
of vibrations to L ) was considered in [7], [J in the framework of Moore's approach 



H and in |12], [15], |28|, |34|, |36| in the framework of the "instantaneous basis" method 
(IBM) described in section 2. The asymptotical solutions to Moore's equation in the 
case suit ^> 1 were obtained in [1(1 EHI, [271, and more general solutions were found in 



jig , |19| , |33| . A detailed study of the problem in the framework of the IBM was given 
in [21] for p = 2 and in for p = 1. The short-time limit eu\t -C 1 for an arbitrary 



integer value of p was considered in [28|. However, in all the cited papers the solutions 



were found under the condition of the strict resonance u w = u v between the mechanical 
and electromagnetic oscillations (excepting recent article ||38|| , where a detuned three- 
dimensional cavity with a nondegenerate spectrum was considered). Evidently, such a 
condition is an idealization. 

The aim of the present paper is to study the case of a nonzero (although small) 
detuning between the frequencies of the mechanical and field modes: 

u) V) =pu x {l + $), |5|<1 (1.1) 

for any integer p = 1,2, . . ., thus generalizing the results of [[H], |22, 28]. It will be shown 



that the photons can be created from vacuum provided the dimensionless detuning 
parameter does not exceed the dimensionless amplitude of the wall vibrations, otherwise 
the total number of photons generated inside the cavity exhibits small oscillations and 
goes periodically to zero. 

The plan of the paper is as follows. In section 2 we give general formulae related 
to the field quantization in a cavity with moving boundaries and derive the simplified 
"reduced equations" in the resonance case. A simple explicit analytical expression for 
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the total energy of the field in all the modes is found in section 3. Section 4 is devoted 
to the "semi-resonance" case p — 1 when the frequency of the wall is close to the 
fundamental frequency of the field. Under this condition new photons are not created, 
but the total energy of all the field modes increases exponentially with time above the 
threshold or oscillates in the case of a large detuning. The generic resonance case of an 
arbitrary p > 2 is analysed in section 5 and the simplifications in the case p = 2 are 
considered in section 6. A brief discussion of the results is given in section 7. Some 
details of calculations are given in the appendix. 



2. Field quantization and reduced equations in the resonance case 



Following the scheme of the field quantization in a cavity with time-dependent boundary 
conditions first proposed by Moore ||, we consider a cavity formed by two infinite ideal 
plates moving in accordance with the prescribed laws 

xieft(t) = u{t), x right (t) = u(t) + L(t) 

where L(t) > is the time dependent length of the cavity. Taking into account 
only the electromagnetic modes whose vector potential is directed along z-axis 
("scalar electrodynamics"), one can write down the field operator in the Heisenberg 
representation A(x, t) at t < (when both the plates were at rest at the positions 
x ieft — and x r i g ht = L ) as (we assume c = h = 1) 



°° \ mrx * 
2 V — = sin — — b n exp (-iu n t) 
In L 



n=l 



h.c. 



(2.1) 



where b n means the usual annihilation photon operator and uo n = ttu/Lq. The choice of 
coefficients in equation ( |2.1| ) corresponds to the standard form of the field Hamiltonian 



H 



±[ L °dx 



57T JO 



dA\ 
dt 



+ 



ox J 



n=l 



(2.2) 



For t > the field operator can be written as 

OO 1 

A{x,t) = -7= \b„ip {n) (x,t) + h.c. 



n=l 

To find the explicit form of functions ^ n '(x, t), n = 1, 2, . . 
that the field operator must satisfy 

i) the wave equation 

(PA _cfA 
dt 2 dx 2 

ii) the boundary conditions 

A(u(t),t) = A{u{t) + L(t),t) = 0, 

iii) the initial condition ( |2.1|) , which is equivalent to 

■?// n) (x, t < 0) = sin exp (—iuj n t) . 



one should take into account 

(2.3) 
(2.4) 
(2.5) 
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Following the approach of Refs. |jl"2 , |i~5| , |I7| we expand the function ifiw^x, t) in a series 
with respect to the instantaneous basis: 

^(x,t > 0) = EQ^Wy^sin ( nk ^ (t)] ) , n = 1,2, . . . (2.6) 

with the initial conditions 

Qfc (0) = 5 fen , Qf '(0) = -iuj n 5 kn , fc,n = l,2, ... 

This way we satisfy automatically both the boundary conditions (|2.4j ) and the initial 
condition (|2]5|). Putting expression (|2.6j ) into the wave equation ( |2.3|) , one can 
arrive after some algebra at an infinite set of coupled differential equations [34, ^ 
(k,n = l,2,...) 

oo oo 

Qi n) + ^tm^ = 2 E + E 9^)Qf ] + O (g%) , (2.7) 

i=i i=i 

where 

w fc (t) = kn/L(t) 

and the time dependent antisymmetric coefficients gkj(t) read (j ^ k) 

, .2/cj [L + ue kj ) . 
to = -9» = (-V ' U 2 - k *)L® > ^ = 1 -(- 1 ) ^ 
For m = (the left wall at rest) the equations like (|2.7|) - ([21f ) were derived in |12|, |I7 |. 

If the wall comes back to its initial position L after some interval of time T, then 
the right-hand side of equation ( |2.7| ) disappears, so at t > T one gets 

Qt\t)=& ) e^ kt + ri { k ) e l ^\ k,n = l,2,... (2.9) 

and ?7j. ri ' ) being some constant coefficients. Consequently, at t > T the initial 
annihilation operators b n cease to be "physical", due to the contribution of the terms 
with "incorrect signs" in the exponentials exp(iu>kt). Introducing a new set of "physical" 
operators a m and a^, which result at t > T in relations such as (|2.1|) and ( |2.2| ), but 
with a m instead of b m , one can easily check that the two sets of operators are related 
by means of the Bogoliubov transformation 

oo 

a m = E (Ka nm + bl(3* nr ^j , m = 1, 2, . . . (2.10) 

n=l 

with the coefficients 

«n m = = (2.H) 

V n V n 

The unitarity of the transformation ( j2.10[ ) implies the following constraints: 



( a nm a fcm ftnmftkm) (^m^ ^m^) (2-12) 

m=l m=l " 

oo oo 

E ( a nm a nj ~ PnmPnj) = E ~ (^"^ £j ^ — ?7m^ ) = ^mj (2.13) 

«=1 n=l " 

OO OO 1 

E go** - p: h a nm ) = e - 0i n M n) - 4 n) U n) ) = o (2.i4) 



n=l n=l 71 
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The mean number of photons in the mth mode equals the average value of the 
operator a^fi™ in the initial state |in) (remember that we use the Heisenberg picture), 
since just this operator has a physical meaning at t > T: 

Af m = (inla^amlin) 

[Pumpkin {pnhk 

n,k 

oo 

+ e ^x^e ] +v^){b%) 



V — lr>)l 2 

/ v I '/m 



\ lilt, I / V / f 



n=l n,ft=l 

+ 2Re E ^M<«- (2-15) 

n,fc=l V^fc 

The first sum in the right-hand sides of the relations above describes the effect of the 
photon creation from vacuum due to the NSCE, while the other sums are different from 
zero only in the case of a nonvacuum initial state of the field. 

To find the coefficients £^ and rj^ one has to solve an infinite set of coupled 
equations ( |2.7|) [k = 1, 2, . . .) with time-dependent coefficients, moreover, each equation 
also contains an infinite number of terms. However, the problem can be essentially 
simplified, if the walls perform small oscillations at the frequency u w close to some 
unperturbed field eigenfrequency: 

L{t) = L (1 + e L sin \pux(l + 5)t]) , u(t) = e u L Q sin [pu^l + S)t + <p] . 

Assuming |ex|, |e u | ~ e 1, it is natural to look for the solutions of equation (|2.7| ) in 
the form similar to fl2.9Q , 

Q%\t) = ^ e ~^+S)t + ^g^M*, (2.16) 

but now we allow the coefficients and r]^ to be slowly varying functions of 
time. The further procedure is well known in the theory of parametrically excited 
systems [41-43]. First we put expression Q2.16Q into equation ( j2.7|) and neglect the 



terms £,77 (having in mind that £,f] ~ e, while £,77 ~ e 2 ), as well as the terms 
proportional to L 2 ~ u 2 ~ e 2 . Multiplying the resulting equation for by the factors 
exp [7x0^(1 + S)t] and exp [— iu k (l + S)t] and performing averaging over fast oscillations 
with the frequencies proportional to Uk (since the functions £, 77 practically do not change 
their values at the time scale of 2ir/uJk) one can verify that only the terms with the 
difference j — k = ±p survive in the right-hand side. Consequently, for even values 
of p the term it in gkj{t) does not make any contribution to the simplified equations 
of motion, thus only the rate of change of the cavity length L/L Q is important in this 
case. On the contrary, if p is an odd number, then the field evolution depends on the 
velocity of the centre of the cavity v c = u + L/2 and does not depend on L alone. These 
interference effects were discussed recently (in the short time limit suit 1) in 



(see also P^j ). We assume hereafter that u = (i.e. that the left wall is at rest), since 
this assumption does not change anything if p is an even number, whereas one should 
simply replace L/L Q by 2v c / L if p is an odd number. 
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The final equations for the coefficients £^ and 77^. contain only three terms with 
simple time independent coefficients in the right-hand sides: 



^i n) = (-l) P [{k + p)& p -(k-p)i 



dr 



(n) 



'-1Y 



{k+p)v, 



(n) 
k+p 



k—p 

(k - p)v ( k n \ 



2i<ykn^\ 



(2.17) 
(2.18) 

(2.19) 

e k "'(0) = Skn, C ; (0) = 0. (2.20) 

Note, however, that uncoupled equations (|2.17|) - (|2.18| ) hold only for k > p. This means 
that they describe the evolution of all the Bogoliubov coefficients only if p = 1. Then 
all the functions T%^(t) are identically equal to zero due to the initial conditions (|2.20|) , 
consequently, no photon can be created from vacuum. If p > 2, we have p — 1 pair of 
coupled equations for the coefficients with lower indices 1 < k < p — 1 



The dimensionless parameters r (a "slow" time) and 7 read (e = el) 
1 

r = -suit, 7 = 6/e. 
The initial conditions are 



dr 



d (n) 

-rVk 
dr fe 



-1? 



(n) 
p—k 



(k + p)^ p - (p - k)r], 
(k + pWkt-(p-k)^ k 



2i"fkr{£ ) . 



(2.21) 
(2.22) 



In this case some functions ^^(t) are not equal to zero at t > 0, thus we have the effect 
of photon creation from the vacuum. 

(n) 

It is convenient to introduce a new set of coefficients p k , whose lower indices run 
over all integers from —00 to 00: 

k > 

k = (2.23) 
k < 



(n) 



0. 



(n) 



Then one can verify that equations (gT7|) - (gT8D and (g^lD - (g^2|) can be combined in 
a single set of equation (A; = ±1, ±2, . . .) 



dr 



Pk 



(n) 



-1^ 



with the initial conditions (n 



Pi n) (0)=4n. 



(k+p)p, 
= 1,2,...) 



(n) 



(A; - p)p 



(») ' 



(2.24) 



(2.25) 



A remarkable feature of the set of equations ( |2.24| ) is that its solutions satisfy exactly 
the unitarity conditions ( |2.12| )-( p7L4"D (although the coefficients £^ and 77^ introduced 
via equation fl2.16p have additional phase factors in comparison with the coefficients 
defined in equation (|2.9|), these phases do not affect the identities concerned), which can 
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be rewritten as 



m=—oo 



(n)* (k) 
m rm 



n8. 



nk 



n, k 



1,2, 



m r (n)* (n) _ (n)* (n) 
Pm Pj P—m P—j 



n=l 

oo 



n 



1,2, 



OO 1 

E 1 W*p { -} 

n L J 



(n)* (n) 
P} P-m 



1,2,... 



(2.26) 
(2.27) 
(2.28) 



n=l 



For example, calculating the derivative I = (d/dr)J2m=-oo m P^m* P^m with the aid of 
equation ( |2.24| ) and its complex conjugated counterpart one can easily verify that 1 = 0. 
Then the value of the right-hand side of ( p. 26 ) is a consequence of the initial conditions 
( p.25| ). The identities fl2.27|) and (|2.28|) can be verified in a similar way, if one uses 
instead of Q2.24| ) the recurrence relations between the coefficients p)$ with the same 
lower index m but with different upper indices derived in section [5]. 

Due to the initial conditions ( |2.25|) the solutions to (|2.24|) satisfy the relation 



pf+2p = if 3 ^ k 

j,fc = 0,l,...,p-l, m = 0,±l,±2, 



(2.29) 



n 



0,1,2,. 



Consequently, the nonzero coefficients p$ form p independent subsets 



ilk — Pj+kp 



(2.30) 



J = 0,1, 



,p-l, g = 0,l,2,.. 
The subset y^' ^ is distinguished, because y^'°^ 



k = 0,±1,±2, 



for k < and the upper index q 
begins at q = 1. This subset is considered in detail in section The generic case is 
studied in section |[ 

3. Total energy and the rate of photon generation 

It is remarkable that to calculate the total energy of the field (normalized by Tioji) 

m 

one does not need explicit expressions of the coefficients p^(r). Calculating the first 
and the second derivatives of £{r) with the aid of the relations ( |2.23| )-( [2T2"B"D one can 
obtain a simple differential equation (see |Appcndix A| ) 



£ = 4p 2 a 2 £ + 4p 2 7 2 £(0) + ^-(p 2 - 1) + 2pVlm(£) 



where 



i -r 



a 



p-i 



Q = 2 E \/ n ( n + P) ( b l b n+ P ) + J2y n (P~ n ) (Kbp-n) 



(3.1) 

(3.2) 
(3.3) 



n=l 



n=l 
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(if p = 1, the last sum in ( |3.3| ) should be replaced by zero). The quantum averaging 
is performed over the initial state of the field (no matter pure or mixed). The initial 
value of the total energy is £(0) = Y^=i n (p\°n) , whereas the initial value of the first 
derivative £(t) reads (see |Appendix A| ) 

£(0) = -paRe(g) (3.4) 

Consequently, the solution to equation ( |3.1| ) can be expressed as 



£(r) = £(0) + 



2 sinh 2 (par) 



£(0) + 



P 



24 2 y ! 



- aRe{g) 



sinh(2par) 
2a ' 



(3.5) 



We see that the total energy increases exponentially at r — > oo, provided 7 < 1. In the 
special case 7 = such asymptotical behaviour of the total energy was obtained also in 
the frameworks of other approaches in [16, 17, |l^, [2^] . Here we have found the explicit 
dependence of the total energy on time in the whole interval < r < 00, as well as 
a nontrivial dependence on the initial state of field, which is contained in the constant 
parameter Q. This parameter is equal to zero for initial Fock or thermal states of the 
field. However, in a generic case Q is different from zero, and it can affect significantly 
the total energy, if £(0) ^> 1. Consider, for example, the case p — 2. If initially the 
first mode (n = 1) was in the coherent state \a) with a = \a\e l( ^, \a\ 3> 1, and all other 

a 2 , so for r ^> 1 and 7 = (exact 
2 e 4r [2 — cos(20)]. The maximal value of the energy 
in this case is three times bigger than the minimal one, depending on the phase 4>. 

According to (pO), the initial stage of the evolution does not depend on the detuning 
parameter 7 for all states which yield lm(Q) = 0, since at r — > one has 

p 2 -l 



modes were not excited, then £(0) = \a\ 2 , Q 
resonance) we have £(r ^> 1) 



£(r) w £(0) - crRe(g) P T + 2 



£(0) + 



24 



(pry 



(3.6) 



Formula ( |3.6|) is exact in the case of 7 = 1. If 7 > 1, then one should replace each 
function sinh(ax)/a in ( p.5|) by its trigonometrical counterpart sin(ax)/a, where 

- 1 (3.7) 

In this case the total energy oscillates in time with the period 7r/(pa), returning to the 
initial value at the end of each period. For a large detuning 7 > 1 the amplitude of 
oscillations decreases as 7 _1 if Re(? 7^ and as 7~ 2 otherwise. For the initial vacuum 
state of field we have 

.2 x 



T 



P 



12a 2 



sinh (par) . 



The total number of photons in all the modes equals M 

00 

J\f(vac) _ 



m,n=l 



n \'lm I 



(3.8) 

j^(vac) + tf(cav)^ where 
(3.9) 
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is the total number of photons generated from vacuum, and the sum 



m,n,fc=l 



m 

Jnk 



'rfc^rfffah) + Re (v^^ (kh))} (3.10) 



describes the influence of the initial state of the field (to obtain the formula (|3.10| ) 
one should take into account the identity (|2.12|) ). Differentiating Q3.9|) and ( p.lOj ) with 
respect to r and performing the summation over m with the help of equations (|2.17| )- 
( |2.22| ) or ( p.24| ) one can obtain the formulae 



dr 
dr 



2aRe - Yj m (P ~ m )p-m{ T )pp-rn{ T ) 



» 



n=l 71 m=l 



(3.11) 



n,k=l 



(n)* (k) (k) (n)=t 

P—m Pp—m i P—mPp—ri 



m=l 



D (Kh) / J (n) (fc) In) (fc) 

2cxRe ^ — =- 2^ m(p - m) p^ip K m L p + p^Pp-r 



n,fc=l 



(3.12) 



m=l 



Consequently, to calculate the total number of photons one has to know the coefficients 
p^ and with the lower indices m = 1,2, . . . ,p — 1. 

4. "Semi-resonance" case (p = 1) 

Let us start calculating the Bogoliubov coefficients with the "semi-resonance" casep = 1. 
It is distinguished, since all the coefficients p£ (t) are equal to zero, and the total number 
of photons is conserved. In this specific case one has to solve the set of equations 
(k,n = l,2,...) 

^ = (* - - (* + + 2ryA*f° (4.1) 

with the initial conditions (0) = <5fc n . To get rid of the infinite number of equations 
we introduce the generating function 



X^\z,T) = Yi { k n \r 



(4.2) 



k=l 



where z is an auxiliary variable. Using the relation kz k = z(dz k /dz) one obtains the 
first-order partial differential equation 



V-i 



2ijz 



dx (n ) 



+ d n) (r 



dr ^ ' dz 

whose solution satisfying the initial condition X^ n \0,z) = z n reads 

z 9{t) - S(t) 



X^{z,r) 



g*(T) 



zS(t) 



+ / &\x)dx 
Jo 



(4.3) 



(4.4) 



where 



S(t) = sinh(ar)/a, g{r) = cosh(ar) +ijS(r) . 



(4.5) 
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e W (T) _ n[-S(r)]^ 
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(4.6) 



[g*(r)} n+1 ' 

Putting this expression into the integral in the right-hand side of equation ( |4.4| j we 
arrive at the final form of the generating function 

z g( T ) - S ( T ) 



g*{r) - zS{t) 



- n 






_ 9*{t) _ 



(4.7) 



which satisfies automatically the necessary boundary condition X^ n \r, 0) = 0. The 
right-hand side of ( |4.7|) can be expanded into the power series of z with the aid of the 
formula ([§4j? vol. 3, section 19.6, equation (16)) 

oo ±m 

(1 - tf- c (l -t + xt)- b =J2 —.{c) m F{-m, b- c; x), 

where F(a, b; c; x) means the Gauss hypergeometric function, and (c)k = T(c + k)/T(c). 
In turn, the function (c) m F(— m, b; c; x) with an integer m is reduced to the Jacobi 
polynomial in accordance with the formula (|44j], vol. 2, section 10.8, equation (16)) 

(c) m F(-m,b;c;x) = m\(-irPt m -^- 1) (2x - 1). 

Consequently, 



(1 - t) b -\l - t + xty h = J2 {-t) m P£- m - c ' c - 1] {2x - 1) 



m=0 



and the coefficient Cm ( T ) rea ds 

£<W( T ) = (_^n-m x n +mp (n-m,-l) ^ _ ^ 



(4- 



(4.9) 



where 



S 



KIT) 



S(r) 



(4.10) 
(4.11) 



s/99 ^1 + S 2 (t) 

K T ) = y/g(T)/g*(T) = \Jl - i 2 k 2 + i7«, | A| 

The form ( fl.9| ) is useful for n > m. To find a convenient formula in the case of n < m 
we introduce the two-dimensional generating function 

oo oo oo 

X(r,z,y) = EEAt'M = T,X {n) (z,r)y n 

yz 



m=l n=l 



n=l 



[g*(r) + yS(r))[g*(r) - g{r)yz + S(r)(y - z)} ' 
The coefficient at z m in (^4.12|) yields another one- dimensional generating function 



(4.12) 



X (ry) = f y^(r) = y ^r)y + S(r)}^ 
Then equation ( P~8| ) results in the expression 

^) = (i _ K 2 ) K m ~«A n+m p^r n,1) (i - 2k 2 ) . 



(4.13) 



(4.14) 
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Note that the functions S(r), cosh(ar) and k(t) are real for any value of 7. For 
7 > 1 it is convenient to use instead of ( |4.5| ) the equivalent expressions in terms of 
the trigonometrical functions: 

S(t) = sin(ar)/a , g(r) = cos(ar) + i-fS(r) . (4-15) 

In the special case 7 = 1 one has S(t) = r and g(r) = 1 + ir. In particular, 

Ui ' 7 ' (i-i r )m+i r n-i [i + r 2 J 1 } 

The knowledge of the two-dimensional generating function enables to verify the 
unitarity condition ( J2.13Q . Consider the product X*(r,zi,yi)X(r,Z2,y2), which is a 
four-variable generating function for the products £ffl*Q ■ Taking y\ = y/uexp(i(p), 
IJ2 = y/u exp(—i(p) and integrating over <p from to 2n one obtains a three- variable 
generating function ^z^z^u™^^*!^. Dividing it by u and integrating the ratio over 
u from to 1 one arrives finally at the relation 

00 1 00 1 

E ^ m 4-d n) *d n) = - In (1 - z\z 2 ) = £ - (zlz 2 ) k , (4.17) 

.-,74 . -1 K 

n,m,l=l k=l 

which is equivalent to the special case of (|2.13|) for 77^ = 0: 

E \&{^f\r) = U mr (4.18) 

Suppose that initially there was a single excited mode labeled with an index n. Due 
to the linearity of the process one may assume that the mean number of photons in this 
mode was u n = 1. Then the mean occupation number of the m-th mode at r > equals 



A/l n) = - = - [(1 - « a )« w| -"P£V 1) (l - ^ 2 )f (4-19) 

where k is given by (gig) . Although formula ( gig) seems asymmetric with respect to 
the indices m and n, actually the relation 

AA") = A^ m) (4.20) 

holds. To prove it we calculate the generating function 

00 

Q(u,v)= E v m u n M™. (4.21) 

m,n=l 

It is related to the function X(z, y) (|4.12| ) as follows 

Q(u, v) = vj- £ dr jf* J* ^X (,/re* ^e*) X* (0^, yfa**) . 
Having performed all the calculations we arrive at the expression 

~, \ 1 + UV — K, 2 (u + V) 

2Q(u, v) = > - 1. (4.22) 

{[1 + UV - K 2 ( U + v)f - AUV(1 - K 2 ) 2 j 

Then ( [4 . 2 0| ) is a consequence of the relation Q{u,v) = Q(v,u). 
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The initial stage of the evolution of N^ij) does not depend on the detuning 
parameter 7, since the principal term of the expansion of (|4.19|) with respect to r yields 



n 



n(n ± 1) . . . (n ± q =f 1) 



2 



However, the further evolution is sensitive to the value of 7. If 7 < 1, then the function 
A/>)(r) has many maxima and minima (especially for large values of m and n), but 
finally it decreases asymptotically as mna 4 / cosh 4 (ar). On the contrary, if 7 > 1, 
then the function M^ij) is periodic with the period ir/a, and it turns into zero for 
t = kn/a, k = 1,2, .. . (excepting the case m = n). The magnitude of the coefficient 
NmH T ) decreases approximately as 7~ 2 l m-ri l for 7 ^> 1. 

In the special case of a cavity filled in with a high-temperature thermal radiation, 
the initial distribution over modes reads v n {T) = T/n, constant T being proportional 
to the temperature. Then Afffi = J2 n ^ni^Afffl. This sum is nothing but T multiplied 
by the coefficient at v m in the Taylor expansion of the function 

" X ^ U r^f 1 ~~ VK 2 (t) 



r 1 c\ii 

Q(v)= / —Q(u,v)=\n 
Jo u 



1-v 
Thus we have 

4 T} = mjVP = T (l - Kt)] 2 ™) . 

We see that the resonance vibrations of the wall cause an effective cooling of the lowest 
electromagnetic modes (provided I7I < 1). The total number of quanta and the total 
energy in this example are formally infinite, due to the equipartition law of the classical 
statistical mechanics. In reality both these quantities are finite, since v n {T) < T/n at 
n —>■ 00 due to the quantum corrections. Other initial conditions in the special case of 
the exact resonance (7 = 0) were considered in |22j . The total energy depends on time 
according to equation (|3.5|) with p = 1. An infinite growth of the energy of a classical 
string whose ends oscillate at the frequency close to U\ in the case of finite amplitude 
and detuning (e ~ 5 ~ 0{\)) was considered in p5|. 



5. Generic resonance case p > 2 

Now we turn to calculating the nonzero Bogoliubov coefficients y^\ T ) ( |2.30| ) in the 
generic case p > 2. One can easily verify that in the distinguished case j = the 
functions y!^' ^ (r) with m > 1 are given by the formulae for £ W (r) found in the preceding 
section, provided one replaces r by apr and 7 by (T7 (remember that a = (— l) p ), 
whereas ym'°K r ) = f° r rn < 0. In the generic case j 7^ it is reasonable to introduce 
a generating function in the form of the Laurent series of an auxiliary variable z 

00 

bM{z,t)= £ y%' j) (r)z m (5.1) 

m=— 00 
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since the lower index of the coefficient yfy 1 ' runs over all integers from — oo to oo. One 
can verify that the function ( |5.1| ) satisfies the homogeneous equation 



Or 



a 



j +pz 



cT 
dz 



zj + 2i>y 

The solution to ( |5.2j ) satisfying the initial condition R^ n ' 3 \z,0) = z n reads 

" zg(pr)+aS(pr)^ n+j/p 



R^\z,t) = z~ j/p 



(5.2) 



(5.3) 



g*(pr) + zaS(pr) 

where the functions S(t) and g{r) were defined in Q4.5|) . The coefficients of the Laurent 
series ( |5.1|) can be calculated with the aid of the Cauchy formula 

' ' dZ -R^\z,r) (5.4) 

in the complex plane in the 



2m Jc z m+1 

where the closed curve C rounds the point z 
counterclockwise direction. Making a scale transformation one can reduce the integral 
( p.4[ ) with the integrand ( |5.3| ) to the integral representation of the Gauss hypergeometric 
function ([[£§, vol 1, section 2.1.3) 

-ir(c) exp(-inb) /•(«+) t b -\l - tf-^ 1 



F(a, 6; c; x) 



-dt, 



(5.5) 



2sin(7r6)r(c- 6)r(6) h (1 - tx) a 

where Re(c — b) > 0, b ^ 1,2,3, .. ., and the integration contour begins at the point 
t = 1 and passes around the point t = in the positive direction. After some algebra 
one can obtain the expression 

{n ,j) = _ r (-m - j/p) r (1 + n + j/p) sin [tt (m + j/p)] 
Vm ' 7rr(l + n-m) 

x (aK) n - m \ m+n+2 J /p F (n + j/p , -m - j/p ; 1 + n - m ; k 2 ) . (5.6) 

We assume hereafter k = k(j>t) and A = A(pr), the functions k(x) and A(x) being 
defined as in ( f4.10| ) and ( |4.11| ). Using the known formula 

T(-z) sin(Trz) = -tt/T(z + 1) (5.7) 

one can eliminate the gamma-function of a negative argument: 



am 



X 



r (1 + n + j/p) (o-K) n - m \ m+n+2 i/P 

r (1 + m + j/p) r (1 + n — m) 
F [n + j /p , —m — j/p; 1 + n — 



m ; k 



(5- 



The form 



gives an explicit expression for the coefficient with < m < n. 



Moreover, it clearly shows the fulfilment of the initial condition y^ ,3 >(r = 0) = 5, 
Transforming the hypergeometric function with the aid of the formula |44], [4 



lim 



F(a,b;c;x) (a) n+1 (b) n+1 x 



n+l 



-F(a + n + 1, b + n + 1; n + 2; x) 



[n 



0,1,2, 



y 



T(c) (n + l)! 

and the identity ( p.7| ) one obtains an equivalent expression 

(n,j) _ 1 ( m +j/P) (-0-K) m ~ n \ m+n+2 i/P 



T (n + j/p) r (1 + m — n) 
x F (m + j/p , —n — j/p ; 1 + m 



(5.9) 
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which gives a convenient form of the coefficient for m > n. Formula (I5.6D with 



vK j) (*>l,n>0): 

(pn+j) _ 
'Ipk-j 



5 j +pm 

negative values of the lower index gives an explicit expression for the nonzero coefficients 

r (k - j/p) r (1 + n + j/p) sin [n (k - j/p)] 

(5.10) 



vrr (1 + n + k) 
x (o-K) n+k \ n - k+2j/p F (n + j/p,k-j/p;l + n + k; k 2 

Note that the expressions ( p.8| )-( |5~TU| ) are valid for j = 0, too. In this case they 
coincide with the formulae obtained in the preceding section. The formulae (|5.8| )-( ^T0| ) 
immediately give the short-time behaviour of the Bogoliubov coefficients at r — > 0: it 
is sufficient to put k « pr, A ~ 1 and to replace the hypergeometric functions by 1. In 
this limit the detuning parameter 7 drops out of the expressions (in the leading terms 
of the Taylor expansions). 

At r — > 00 we have the following asymptotics of the functions k(pt) and X(pr) (if 
7<1) 

1 



1 



A — ► a + ry, r 



00. 



Then equation (|5.6| ) together with the known asymptotics of the hypergeometric function 
F(a, b; a + b + 1; 1 - x) at x < 1 [g [§§ 



F(a, &; a + 6 + 1; 1 — x) 



lead to the asymptotical expression for the Bogoliubov coefficients 

sin[7r(m + j/p)] 



r > 1) 



ir(m + j/p) 



-(a + ry) 



i + o(^.„s 



(5.12) 



For 7 < 1 the correction has an order mnr exp(—2apr) , while for 7 = 1 it has an order 
mnln(r)/r 2 . 

One can verify that the generating function ( |5.3| ) satisfies the recurrence relation 

= (j + q p) | a [fl(9-^) - i^+U)] + 2z 7j R (,?j) } (5.13) 

Its immediate consequence is an analogous relation for the Bogoliubov coefficients with 
the same lower indices: 



7 r m 

get 



{> 



n < a 



Jn-p) 
rm 



n ( n +P) 
rm 



(5.14) 



Equation ( |5 . 1 4] ) is valid for n > p (when q > 1 and j > 1 in ( |5.13| )), since the coefficients 
are not defined when n < 0. However, using the chain of identities 



1 r 



-j/p 



z L 



#(o,p-i)(l/V) 



5 + 

1 



j/p-i 



E vl 



(0,p-j)* 









[ <? + 5/2 J 




00 

E vS£ 
fc=— 00 



)* ~k 
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one can obtain the first p — 1 recurrence relations 



d_ 

dr 1 



» 



n < a 



(p—n)* 



rra 



n 



l,2,...,p-l. (5.15) 



To treat the special case n — p (it corresponds to the distinguished subset with j = 0) 
one should take into account that Bj°'°\z) = 1, which means formally that p^jjj = <5 m o- 
So the last recurrence relation reads 

±pM= p {-*pto) + 2y rp <p)} t m >l 

(remember that p'jjj = for m < 0). Now one can verify that the unitarity conditions 
( p.27| )-( p~2^ ) are the consequencies of the equations ( |5.14|) and (|5.15| ). 

Differentiating the "vacuum" part of sum ( |2.15|) with respect to r and performing 
the summation over the upper index n with the aid of (|5.14j )- (|5.15| ) (remembering that 
the coefficients p^jj/j are different from zero provided the difference n — m is a multiple 
of p) one can obtain the formula for the photon generation rate from vacuum in each 
mode (0 < j < p - 1, q = 0, 1, 2, . . .) 
d 



Kflyac) 

dr 3+pq 



2W1 



-2a(j+pq)Re[^ pq r ] ! J p -J q ) 

— sm(irj/p)T(q + j/p)T(l + q + j/p)T(2 - j/p) 2qA 
TrT{j/p)T{q+l)T{q + 2) 



x 



F (q + j/p, -j/p; 1 + q ; k 2 ) F (q + j/p , 1 - j/p ; 2 + q ; k 2 ) (5.16) 

We see that there is no photon creation in the modes with numbers p, 2p, .... At r 1 
we have N^/jjjq ~ r 2q+1 . In the long-time limit the photon generation rate tends to the 
constant value (if 7 < 1) 



d ir{vac) 

dr 3+pq 



2ap 2 sin 2 (it j/p) 



apr 3> 1 



(5.17) 



ir 2 (j +pq) 

For g > 1 and for a fixed value of k one can simplify the right-hand side of (|5.16|) using 
Stirling's formula for the Gamma-functions and the easily verified asymptotical formula 



In this case 



F(a, b; c; z) ~ (1 — az/c)~ 



d_ ir(vac) 

dr J+pq 



2pJl 



7 2 k 2 



a,c>l. 

sin(7rj/p)r(2 - j/p)K 2q+l 
nTij/p^-iM (I - K 2 ) 1 - 2j/p ' 



g>l. (5.U 



In particular, if q ^> S 2 (pr)\ ^> 1, then 



2 P a M^MT^mif/qr^ exp ( _ g/ ^ 2 



(5.19) 



7rT(j/p)S 2 

Comparing (|5.17| ) and (|5.19| ) one can conclude that the number of the effectively excited 
modes (i.e. the modes with a time independent photon generation rate) increases in time 
exponentially, approximately as S 2 (r)/ In S(r). 
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Differentiating equation ( |3.11| ) once again over r one can perform the summation 
over the upper index n with the aid of equations (|5.14|) - (p.l5|) to obtain a closed 
expression for the second derivative of the total number of "vacuum" photons 



J2 p-i 

Kf(vac) = 2Re y / 
m=l 



t(m)Ap-m) 
Sm Sp— m 



i „(p-"i)*„( m )* 
^ '/m ' Ip—m 



p-1 

2 ^ m(p — m) < m(p — m) 

m=l 



+ (l - 2 7 V) F I — , ; 1 ; « 




In the short-time limit one obtains 
M {vac) = \p{p 2 - 1), 



\apr\ <C 1 



In the long-time limit the formulae ( |5.7| ), ( |5.11| ) and Z)m=i sm (vrm/p) 
another simple expression (provided p > 2) 

tf(vac) = 2aV/7r 2 , apr > 1, a > 



(5.20) 



(5.21) 
p/2 lead to 

(5.22) 



Consequently, the total number of photons created from vacuum due to NSCE increases 
in time quadratically both in the short-time and in the long-time limits (although with 
different coefficients). 



It is interesting to compare formula (|5.22|) with the total rate of change of the 
number of "cavity" photons due to nonvacuum initial conditions. Using equation ( |3.12| ) 
and replacing the coefficients pffl by their asymptotical values (|5.12|) one can obtain the 
expression 



d ^f(cav) Aap 



dr 



7T 



2 P-1 
m=l 



sin 



(nm/p) 

n,k=0 



a 



n+k 



m + pn)(m + pk) 



^m+prfim+pk] 



) (a + 2 7 )^ n - ^Re [(b m+pn b m+pk ) (a + * 7 ) fc+ " +1 ] } (5.23) 



which holds provided apr ^> 1 and a > 0. For the physical initial states the sum in 



the right-hand side of ( 5.23 ) is finite. This is obvious if a finite number of modes was 
excited initially. But even if the cavity was initially in a high-temperature thermal 
state, so that (b\fik) = o~nkT/n, (b n bk) = 0, the sum over n,k yields a finite value 
TJ2^=o {iTi + pn)~ 2 . Consequently, the total number of "nonvacuum" photons increases 
in time linearly at apr ^> 1, whereas the total number of quanta generated from vacuum 
increases quadratically in the long time limit. At the same time, the total "vacuum" and 
"nonvacuum" energies increase exponentially if 7 < 1 (see section 3). The origin of the 
difference in the behaviours of the total energy and the total number of photons becomes 
clear, if one looks at the asymptotical formulae ( |5.17| )-( |5.19| ). They show that the rate of 
photon generation in the mth completely excited mode decreases approximately as 1/m 
(excepting the modes whose numbers are multiples of p), so the stationary rate of the 
energy generation asymptotically almost does not depend on m. In turn, the number 



Resonance photon generation in a vibrating cavity 



17 



of the effectively excited modes increases in time exponentially. These two factors lead 
to the exponential growth of the total energy (see also [E7J in the special case 7 = 0). 



6. The "principal resonance" (p = 2) 

Some formulae obtained in the preceding section can be simplified in the special case 
p = 2. In this case there are two subsets of nonzero Bogoliubov coefficients. The first 
one consists of the coefficients with even upper and lower indices t^f* which are reduced 
to the coefficients ^ of the "semi-resonance" case (since = 0, this subset does 
not contribute to the generation of new photons). The second subset is formed by the 
"odd" coefficients which can be written as [k = k(2t)] 

r(m L Q /0\ ty-Ti— m \ m+n+1 
,_„ , , [n + 6/1) k A 



e 



2 "' +1 r(m + 3/2)r(l +n-m) 

x F (n + 1/2 , -m -1/2; l + n-m; k 2 

(2n+i) _ (-l) m ~ n r (m + 1/2) K ™-«A m+n+1 



n > m 



(6.1) 



2m+l 



'hk+1 



r(n + l/2)r(l + m-n) 
F (m + 1/2 , -n - 1/2 ; 1 + m-n; k 2 



m > n 



(6.2) 



n+fc+l^n— k 



(-l) fc -T (k + 1/2) T (n + 3/2) «» 
?rr (2 + n + k) 

F (n + 1/2 , k + 1/2 ; 2 + n + fe ; k 2 ) . 



(6.3) 



All the "odd" coefficients can be expressed in terms of the complete elliptic integrals 

tt/2 



47 



K k) 



da 



yl — k 2 sin 2 a 



E K 



vr/2 



dav 1 — k 2 sin a . 



In particular, 



where 



(i) 



7T 



-A(«)E(«), n\ 



(i) 



7VK 



R 2 K(k)-E(k) 



1 + S 2 {2t) 



-1/2 



(6.4) 



(6.5) 



However, the analogous expressions for the coefficients £W and r/^^ 1 with m, n > 1 
appear rather cumbersome (they can be written as linear combinations of the functions 
E(k) and K(k) multiplied by some rational functions of K and R), so we do not bring 
them here. 

The photon generation rate from vacuum in the principal cavity mode (m = 1) 
reads 

dMi vac) 



dr 



-2Re 



8Vi -i 2 n 2 



it 2 k 



E (k) 



E(k) - k 2 K(K) 



(6.6) 
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The total number of photons in the first mode can be obtained by integrating equation 
(|6.6D. Taking into account the relation 

T - 7 2 K 2 dr = dn/k 2 (6.7) 

(6.8) 
(6.9) 



and the differentiation rules for the complete elliptic integrals 
dK(«) E(k) K(k) dE(«) E(k)-K(k) 



dn kk 2 k 
one can verify the following result: 

7T Z 



dK 



A/i (mc) («) 



2E(«) - k 2 K(k 



Making the transformation |44], |46 
1 



K 



1 + k 



K K 



E 



ft 



ft 



1 

2" 

E(«) + kK(k) 



one can rewrite formulae (|6.4j) and (|6.9|) in the form given in [21j for 7 = 0. Using the 



asymptotical expansions of the elliptic integrals at k — » 1 |48 



E(«) 



In - + - ( In - - 1 I k 2 + 

ft 4 V ft 

1 /, 4 1' 



In 



ft 



R 2 + -. 



one can obtain the formula 

A/? ac) (T> 1) 



8a 4/2 
— r + — In - 

7T Z 7T Va 



i + O (re 



-4ar 



a > 0. 



(6.10) 



In the special case of 7 = 1 one can obtain the expansion 



(vac) 



T > 1) 



4 , 12 1 
— In r + — In 2 



O 



7T^ 



It 



_-2 



T 



If 7 > 1, the number of photons in the principal mode oscillates with the period 7r/(2a) 
For 7 ^> 1 one can write ft rs sin(2ar)/a, i.e. |ft| <C 1. In this case 

k 2 sin 2 (2ar) 

The second derivative of the total number of "vacuum" photons can be written as 



<c 1. 



dr 2 



Re 



(i) 



+ 



*7i 



(i) 



it 2 k 2 



(6.11) 



ft 4 K 2 (ft) - 2ft 2 K(ft)E(ft) + (l + k 2 - 2 7 2 ft 4 ) E 2 (ft) 
In the limiting cases this formula yields 

AA mc )(r< 1) ^r 2 

Af (mc) (r> 1) = 8a 2 r 2 /7r 2 + C(r), a > 0. 

If 7 3> 1, then |«| <C 1, but 7 2 ft 2 « sin 2 (2ar) ~ In this case the Taylor 

expansion of the expression ( p. 11 ) yields J\f^ vac ^ = 2 cos(45t) + 0(-~f~ 2 ). Integrating this 
equation with account of the initial conditions J\f( vac \0) = J\f( vac \0) = one obtains 
M {vac) K tfiv°c) ~ sin 2 (2&r)/(4S 2 ). 
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7. Discussion 

Let us discuss briefly the main results of the paper. We have solved the problem of 
the photon generation due to the nonstationary Casimir effect in an ideal Fabry-Perot 
cavity with an equidistant spectrum, if the cavity walls perform small (quasi)resonance 

oscillations at the frequency u w = p(ttc/ Lq)(1 + 5), for any integer value of p — 1, 2, 

Namely, we have found explicit analytical expressions for the Bogoliubov coefficients, 
the rate of photon production in each mode and the total energy in the case of an 
arbitrary (although small compared with u w ) detuning. These expressions are exact 
consequences of the reduced equations (|2.24j) or ( |5.14|) -( ^T5|) . One should remember, 



however, that the reduced equations arise after averaging the exact equation (|2]7|) over 
fast oscillations and neglecting the second-order terms with respect to small parameters 
e and 5. Consequently, the "true" functions N(t), £(t), etc. could differ from those given 
above in terms proportional to e 2 . But such a difference seems quite insignificant under 



the realistic conditions. As was shown in [18, 21], it is hardly possible to obtain the 



value of the dimensionless amplitude of the resonance wall vibrations e exceeding 10~ 8 
in a laboratory. This means that the relative difference between the "true" magnitude 
of the photon generation rate (for example) and that given in section 5 could be of the 
order of 10 -8 (or less) for t < t c ~ (c^e 2 ) -1 . For uj\ ~ 10 10 s _1 the characteristic time t c 
has an order of months or years, and even for the optical frequences it has an order of 
seconds (although it is unclear how to cause the wall to vibrate at an optical frequency 
with a sufficiently big amplitude). Another argument in favour of the solutions obtained 
is that these solutions satisfy exactly the Bogoliubov transformation unitarity conditions 

O)-(Hl)- 

Note that the rate of photon generation from vacuum in some mode is proportional 
to p 2 e (if 7 = 0), and the total generation rate is proportional to p 3 e 2 . Actually, 
the dimensionless amplitude of the wall oscillations e is inversly proportional to the 



frequency, since it is determined by the maximal possible stresses inside the wall |18|, [21 
Thus we see that increasing the resonance frequency one could achieve, in principle, some 
amplification of the number of photons proportional to p. 

It was shown in the previous studies [10,14-23] that the photon production from 
vacuum due to the NSCE could be observed under the condition of the strict parametric 
resonance. Here it is demonstrated explicitly that the photons cannot be produced if 
the detuning S exceeds the dimensionless amplitude e. This result confirms once again 



the statement made in [^1|] that the NSCE could be observed only in the resonance 
regime, ruling out the nonresonance laws of motion of the wall. The requirements to 
a possible experiment turn out rather hard (for example, for the principal frequency 
about 10 GHz the detuning should not exceed 100 Hz for the time of the order of at 
least 0.01 s), but they do not seem to be absolutely unrealizable. 

Another source of troubles is connected with a nonideality of real cavities. Until 
now there were only few attempts to take into account different losses in the cavities with 
moving boundaries p3[ |3l| , and this problem is still a challenge for theoreticians. 
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Using equations ( |2.15| ) and ( |2.23| ) one can express 

OO 

n 



S = J2 ~S {n) + f] &Mui nk) + Re E ^T^-t-o 



where 



OO 



m 



1=1 



n,k=l 



in) 
P-m 



the total energy in all the modes as 
{b n bk) TT {nk) 

n,k=l 



OO OO 

m= — oo m.= — i 



2 n (") n ( fc ) 



-in 



(A.l) 
(A.2) 



m=— oo 



(to write as a sum from — oo to oo one should take into account that the summand 
in the last sum of (|2.15 ) is symmetrical with respect to n and k). Differentiating Jj[ nk ^ 
with respect to r and taking into account the equations fl2.24|) one can obtain after a 
simple algebra the expression 



m(m 



J oo 
aT m=-oo 

Differentiating the above expression once more one obtains 
E#*> = Ap 2 u[ nk) + 2z 7 p 2 (-l)P X S" fc) , 

where 

oo 

(nk) _ \- / x [ (fc) (n)* _ (n)* (fc) 

CI — / , m \ m ^ P) [Pm Pm+p Pm Pm+p 
m,= — oo 



rm rm+p 1 rm Hm+p 



(A.3) 



Differentiating xi one can verify that 

= 2i 1 {-iyu[ nk) . 

Consequently, 

x [" fc ) = 2* 7 (-l) p f/ 1 (nfc) + const, 

where the additive constant is determined by the initial conditions. Finally we arrive 
at the equation 

u[ nk) = V(i - 7 2 )f/i (nfc) + V 7 2 ^i (nfe) (o) + 2i 1 p 2 {-iy x ( i k \v) 



with 



ll[ nk \o) = n 2 5 nk , Xi" fe) (0) = nk [5 kyn - p - 5 n>k - P } ■ 



Using the same scheme one can obtain analogous relations for the coefficient U% . 



J oo 

^Uf k) =p(-IT E rnp^{(m + P ) P ^L p -( P -m) Pp k 2 



m=—oo 



dr 

= Ap 2 ui nk) - 2i 1 p 2 (-lf X { 2 k \ 



(A.4) 
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X2 k) = m Pm [i m + P)p-m- P + (P - fn)pp-m 

m=—co 

X { t ] = -2* 7 (-im (nfc) , 

Uf k) = 4p 2 (l - 7 2 )f/ 2 (nfe) - 2z 7 p 2 (-l)^ ) (0), 
X2 (°) = nk5 k)P - n . 
The calculation of the vacuum contribution to the total energy 



"fin 



n=l 



is more involved, since the summation in ( |A.1|) is performed now not from — oo to oo, 
but over the coefficients pffl with negative indices m only. Differentiating the sum ( A.l|) 
with respect to r and using equations ( |2.24|) we obtain 

oo 

S(n) = 2 (_l)PRe X: m 2 pLl [(m + p)p^ p + (p - m)p^] . (A.5) 

m=l 

Differentiating the expression ( |A.5|) once again one can obtain after some algebra the 
equation 

OO I 



n=l 



m=l 



+ 2Re V] m(p — m) 2 [mF m + (m + p)G r , 



(A.6) 



where 



$W = Im ™?P-m 



(p - m)p p n l* m 



(m + p)p^t- P 



m=l 



oo 1 

f = V - L (n) V n) - 

— i ' i 



(n)* (n) 
P—m P—m 



n=l 
oo 



OO 1 

™ / y ^"m+pfVt— p Pp—mr—m—p 



n=l 



Differentiating the function over r and using again equations (|2.24j) one can verify 
that the derivative d^/dr of the combination \P = Y^=\ - $ (n) + p(-l) p 7>S' (n) can be 
written in the form analogous to the last sum (from 1 to p) of equation (|A.6|) , but the 
symbol Re should be replaced by Im. Since F m = 1/m due to the identity (|2.27|) and 
G m = due to (|2.28|) , we have d^/dr = 0. Taking into account the initial conditions 
$( n )(0) = S( n '(0) = one obtains ^(t) = 0. Combining all the terms giving the second 
derivative of £ one can arrive finally at equation ( |3.1|) , where the term \p 1 {j? — 1) is 
the value of the sum 2 YZn=\ m{p — m) 2 . 

The initial value of the first derivative £{t) is determined by the right-hand sides 
of equations (|A.3|) , QA.4Q and (|A.5|) taken at r = 0, when pffl = 5 mn : 

£ (0) = -2pa J2\ n ( n + P)^ e i b l b n+ P ) -po-J2\ n (P~ n)Re{b n b p _ n ) . 

n=l n=l 

Comparing this formula with (|3.3| ) we arrive at equation 
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